Theory of the spin relaxation of conduction electrons in silicon 
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A realistic pseudopotential model is introduced to investigate the phonon-induced spin relaxation 
of conduction electrons in bulk silicon. We find a surprisingly subtle interference of the Elliott 
and Yafet processes affecting the spin relaxation over a wide temperature range, suppressing the 
significance of the intravalley spin-flip scattering, previously considered dominant, above roughly 
120 K. The calculated spin relaxation times Ti agree with the spin resonance and spin injection 
data, following a T -3 temperature dependence. The valley anisotropy of Ti and the spin relaxation 
rates for hot electrons are predicted. 

PACS numbers: 72.25.Rb, 72.25.Dc, 76.30.Pk 



Silicon is the core material for the information tech- 
nology. Yet we know surprisingly little about the spin 
relaxation processes of its conduction electrons P, Q • 
With the pioneering demonstration of the spin injection 
into silicon by Appelbaum et al. 0, 0] and related ex- 
perimental breakthroughs 0, 0] and steps towards sili- 
con s pin tronics Q, as well as with theoretical analyses 
@, SUB], we have gained valuable new insight into the 
spin transport and relaxation in this material. 

A systematic investigation of the conduction electron 
spin relaxation time Ti in silicon was conducted earlier 
by Lepine using electron spin resonance (ESR) (see 
the summary in Ref. 0). As the ESR does not discrim- 
inate between conduction electrons and electrons bound 
on donors, the representative data are limited to temper- 
atures above T w 150 K at which most electrons are in 
the conduction band for the investigated donor densities 
111] , Spin injection, on the other hand, looks at conduc- 
tion electrons only. Appelbaum et al. d, 0] extracted 
useful data below 150 K, filling the gap. These mea- 
surements were using samples with nondegenerate elec- 
tron densities. The conduction electron spin relaxation 
and spin transport properties were also investigated in 
Si/SiGe quantum wells [1(1 [Hj]. There, however, the spin 
coherence is due not to the bulk-derived properties, but 
rather due to the appearance of the structure-inversion 
anisotropy spin-orbit fields. 

While it is generally believed that the spin relax- 
ation in silicon is caused by the Elliott- Yafet mechanism 
0,1,0, ED, El (unlike in III-V semiconductors, in which 
the most important mechanism is the D'yakonov-Perel' 
one there is yet no systematic theoretical study of 
it. There are two processes involved: that of Elliott 
and Yafet. In the Elliott processes 15[ the spin- flip is 
due to the admixture of the Pauli up and down spins 
in the Bloch state, caused by spin-orbit coupling. The 
electron-phonon matrix element couples only equal Pauli 
spins. In the Yafet processes [l3| spin flips are due to 
the phonon-modulated spin-orbit coupling so that the 
electron-phonon coupling alone couples opposite spins. 



The two processes interfere destructively at low phonon 
momenta affecting Ti typically at very low temperatures 
[l3j |. Yafet gave qualitative estimates for T\ in silicon 
assuming intravalley electron-acoustic phonon scatter- 
ing, finding that T\ ~ T -5 / 2 . This temperature de- 
pendence has been widely used to fit experimental data 
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Here we perform comprehensive theoretical investiga- 
tion of Ti in bulk silicon within the Elliott- Yafet mech- 
anism. We introduce a pseudopotential model that re- 
produces the known spin-orbit splittings of the relevant 
electronic states. The model, together with a realistic 
phonon structure taken from the adiabatic bond charge 
model 17], allows us to calculate the spin mixing prob- 
abilities and the electron-phonon-induced spin flips for 
both the Elliott and Yafet processes. We show that the 
interference between the two processes affects T\ over 
a remarkably wide temperature range. Both the intra 
and intervalley spin-flip scatterings are important [opti- 
cal (OP) phonons are less relevant than acoustic (AC) 
ones] making Yafet 's prediction invalid. Our calculated 
Ti(T), which is in quantitative agreement with experi- 
ment, is well described by the T\ ~ T~ 3 dependence. 
We further predict the valley anisotropy of T\ and give 
the spin relaxation rates for hot electrons. 

The presence of space inversion symmetry in bulk sil- 
icon allows to write the Bloch states as combinations of 
the Pauli spinors @, El , 
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Here k is the lattice momentum confined to the first Bril- 
louin Zone, n is the band index, 'fr(-IJ-) is the effective spin 
index, g denote the reciprocal lattice vectors, and |k) 
stand for the plane waves. The two states above are 
degenerate, and can so be chosen as the spin "up" and 
"down" states to satisfy (k, n -ff \a z \k, n JJ,) = 0. The mix- 
ing of the Pauli spins in the spin "up" ( "down" ) state is 



2 



characterized by the mixing probability |&k,n| 2 , which is 
key to understand spin relaxation. 

We follow the scheme of Rcf. and build a pseu- 
dopotential model incorporating spin-orbit coupling to 
obtain the electronic states and electron-phonon cou- 
pling needed to calculate T\. The pseudopotential of 
each atom is the sum of the scalar (y n ) and spin-orbit 
(v so ) parts, v(r) = v n (r) + v so (r); the form factors 
are u(ki,k 2 ) = J(d 3 r/a 3 )e-^ kl " k2 > r w(r), with the lat- 
tice constant a = 5.431 A. The scalar part v n (r) is 
taken from Ref. [l||, reproducing well the silicon bands. 



21) v so (r) 



We introduce the spin-orbit part as [2 
X9(r— r c )Iu-&Vi, where L is the angular momentum oper- 
ator, Vi is the projector on the orbital momentum state I, 
r c = 2tb is the effective radius for the spin-orbit coupling 
(tb is the Bohr radius), and A is a spin-orbit parameter 
to be fitted to the valence band spin-orbit splitting A_L . 
The fitted form factor of v so can be approximated by 21 1 
u so (ki,k 2 ) = — iA so (a/27r) 2 ki x k 2 ■ <x, with the effec- 
tive spin-orbit interaction A so = (167r 3 /45)(r c /a) 5 A = 
9.475 x 10 -4 eV. The actual spin-orbit pseudopotential 
Hamiltonian matrix element comprises the structure and 
form factors: (ki|i? so |k 2 ) = cos[(ki - k 2 ) • r]w so (ki, k 2 ), 
where r = (a/8)(l, 1, 1). The pseudopotential Hamilto- 
nian is diagonalized on a basis of 387 plane- waves, large 
enough to converge the spin mixing probabilities around 
the conduction band edge. The sums over momenta are 
performed using the tetrahedron method. 
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FIG. 1: Calculated electronic band structure of silicon. The 
single group notation 22] is used in the main plot, while the 
double group notation [231 ] appears on the zoom plots. 

The calculated electronic band structure is shown 
in Fig. [TJ Table Q] displays selected band-structure 
properties, calculated and measured. The agreement 
with known data is very satisfactory, justifying our 
pseudopotential for exploring spin-orbit effects in sil- 
icon. In the following we denote the wave vectors 
of the six conduction-band valleys as K^, with i — 
X(X)/Y(Y)/Z(Z) standing for the corresponding valley 
orientation. 

Let us first see what can we learn about b 2 from sym- 
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Exp. [24] 
Cal. 


3.4 1.25 1.17 
3.46 1.33 1.19 


44 30-40 
44 35 


0.9163 0.1905 
0.915 0.204 


0.846 



TABLE I: Comparison between measured and calculated 
band-structure characteristics. Displayed is the energy Er 15 
of the direct gap at Fi5, the energy Ex ± of the Xi (Xs) point, 
and the indirect band gap E g . The spin-orbit split off energy 
of the top of the valence band is A25, and that of the con- 
duction band Fis point is A15 (the stated range of values are 
calculated 24]). Further shown are the longitudinal (mf 1 ) 
and transverse (m^) effective masses (in the units of the free 
electron mass) at the conduction band minimum at fenm along 
the A lines FX. 



metry group arguments. Consider the Z valley and take 
the spin-orbit interaction H so cx L x a x + L y a y + L z a z as 
a perturbation. Group theory shows that the conduction 
band (Ai, see Fig. [l| couples only to the the valence 
band (A5) at the band edge, so that the spin mixing is of 



the order of |6 cfc | 2 » (A so /E g ) 2 - 1CT 6 [I3j,[l5j. In addi- 
tion to the magnitude, we can also learn about the valley 
anisotropy of b 2 . We can describe the symmetry charac- 
ter of the conduction band orbital states as \Z), and that 
of the degenerate valence states as |X)/|Y) [23j, in the 
usual kp theory sense. The only nonvanishing matrix el- 
ements of the orbital momentum between these states are 
(X\L y \Z) = -(Y\L X \Z), by symmetry. Therefore, the ef- 
fective spin-orbit interaction involves only the a x and a y 
terms with equal contribution. If the spin is quantized 
along z, both a x and a y terms contribute equally to the 
spin mixing b 2 . However, only one of them, i.e., a y /a x 
term, contributes to the spin mixing for the x/y quan- 
tized spin. Since the spin along x/y for the Z valley is 
equivalent to the spin along z for the X/Y valley under 
a ^-rotation around the y/x axis, we can conclude that 
the average spin mixing is anisotropic with respect to the 
valley orientation: (& 2 ) £ ,z ~ ^(b 2 ) e ,x/Y- This is indeed 
found from numerics, as shown in Fig. [21 See EPAPS 
Document No. for a more intuitive explanation of the 
anisotropy. 

Figure [2] also shows a peak around 150 meV, which 
is roughly the energy of X\. This peak indicates a 
large spin mixing on the X plane. In the X plane, 
there is anticrossing along the X-W direction due to 
spin-orbit coupling of degenerate bands Ai and A 2 , 
which results in large spin mixing (similar to spin hot 
spots 20]). The anticrossing is absent at the X-point, 
but grows along the X-W direction to the order of | 
quickly. See EPAPS Document No. for more details 
on spin hot spots in silicon. In the inset of Fig. [2 
we also plot the temperature dependence of the aver- 
age spin mixing {b 2 ) T ,i = Ekei \ b ck\ 2 f (^k) / Ekei /( e k), 
where /(ek) = Ce £k ^ feflT ^ is the properly normalized 
(C) Maxwell-Boltzmann distribution. The anisotropy, 
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FIG. 2: Calculated energy dependence of the spin mixing 
probability (6 2 ) e in the X/Y (dashed) and Z (solid) valleys. 
The dot-dashed curve is 2{b 2 ) € ^x- The inset shows the tem- 
perature dependent (6 2 )t- 



)t,z ~ 2(b 2 )T,x/Yi remains. 
The spin-flip electron-phonon matrix elements for the 
conduction states |ki;cf|*) and |k2 ; c J|) are [lj| 



M A (k 1) k 2 ) 




q glg2 



x( fl ^ c(gl) l ,(k 1 + gl ,k 2 + g2) h b ^M) 



v ^k 1)C (g2) 



(g 2 ) 



,(3) 



where q = ki — k 2 is the phonon wave vector, r a = ±t 
are the position vectors of the two basis atoms, Ak = 
ki + gi — k 2 — g 2 , p is the silicon density, and h£l^ and 
£q a are the phonon energy and the polarization vector 
for the evth atom, obtained here from the adiabatic bond- 
charge model [17] . For the Elliott processes v = v n , while 
for Yafet processes v = v so in Eq. ([3]). 

The spin relaxation rate is given by T± (T) = 

Ei,i r iT«> with 



T iM T ) = H J deT*(e,T)f(e). 



(4) 



Here T^- is the total scattering rate for the electrons of 
energy e making A-phonon assisted spin-flip transitions 
from the ith to the jth valley. For non-degenerate elec- 
tron densities, 



x (n qA + ^±i)(5(£ ck2 -e±^ qA ), (5) 

where fiq,x is the phonon distribution. Considering the 
low-energy intravalley processes with AC phonons Yafet 



found the dependence of T| ntra - T~ 5 / 2 [13]. That pre- 
diction was based on Yafet 's observation that the time 
reversal and space inversion symmetry inhibits the low 
momentum scattering, expressed in the quadratic depen- 
dence, M A (ki,k 2 ) oc |q| 2 /wfi^, valid for low q. 
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FIG. 3: Spin relaxation time in silicon as a function of tem- 
perature. The solid curve is the calculation, the symbols are 
the spin injection (Appelbaum) I3l, 12511 and the electron spin 
resonance (Lepine and Lancaster) [ill, flrjl] data (see Ref. Q). 
The dashed-dotted curve is J n 1 lntra of the intravalley scattering 
only. The inset shows the contour plot of the spin relaxation 
rate of hot electrons, as a function of the electron energy 
e and lattice temperature T. 

Our main result is in Fig. [3]which shows the calculated 
Ti(T). The agreement with experiments @, i, [n|, [l6| 
is very good. Clearly, the intravalley spin relaxation 
juntra a i one j s n0 ^ sufficient to explain the experiment. 
While there is no single scattering type governing the 
whole temperature range, precluding a simple theoret- 
ical prediction for Ti(T), a fit to the numerical data 
gives T\ ~ T~ 3 in the investigated temperature region. 
Considering intravalley processes only, Yafet 's prediction 



T 5 / 2 works at low T. Due to the strong spec- 



tral dependence (peak) of (& 2 ) e more energetic electrons 
at higher temperatures strongly bias the average spin re- 
laxation, making y i lntra (T') decay faster than predicted. 
This behavior is also reflected in the energy-resolved spin 
relaxation rate, shown in the inset to Fig. [31 The strong 
increase of the rate with increasing energy is due to the 
increase in the scattering phase space and b 2 . 

Figure |4] gives a comprehensive analysis of the calcu- 
lated T\. The comparison of the AC and OP phonon 
contributions as well as the intra and intervalley contri- 
butions to spin relaxation is shown in Fig. BJa). We 
find that the spin relaxation is dominated by intraval- 
ley electron-AC phonon scattering at low temperatures 
(T < 120 K), and by intervalley electron-AC phonon 
scattering at high temperatures (T > 120 K). Figure 
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FIG. 4: Analysis of the spin relaxation in silicon, (a) Spin 
relaxation time induced by acoustic (green solid) and opti- 
cal (green dashed) phonons, as well as by intervalley (dash- 
dotted) and intravalley (solid) scattering, (b) Elliott (green) 
and Yafet (red) processes compared to the total Ti (blue). 
The estimate based on the knowledge of the momentum re- 
laxation time r and b 2 is also indicated (dashed blue) . (c) El- 
liott and Yafet processes with and without intravalley acoustic 
phonons. The blue curve is the total T\. (d) Valley-anisotropy 
ratios Tix/Tz (solid); intravalley (dashed) and intervalley 
(dot-dashed) scattering are distinguished. 



HJb) resolves the individual contributions of the Elliott 
and Yafet processes. Remarkably, they would individ- 
ually lead to Ti orders of magnitude below the actual 
spin lifetime, over the whole temperature range! The de- 
structive interference of the two processes is a very subtle 
effect (also in terms of numerics). The fact that the El- 
liott and Yafet processes interfere destructively was pre- 
dicted already by Yafet on symmetry grounds [13j ]. See 
EPAPS Document No. for more details. The same fig- 
ure also brings the estimated spin relaxation time using 
the conventional formula 0,13 Ti = T/4(b 2 ) T , where 
r is the momentum relaxation time calculated from our 
pseudopotential model. This estimate fails in silicon es- 
pecially at high temperatures. The subtle nature of the 
interference between the Elliott and Yafet processes is 
displayed in Fig. IDJc). Each process is dominated by the 
the intravalley electron-AC phonon scattering over all 
temperatures. The interference drastically reduces the 
significance of this scattering. One consequence of our 
calculation would be the decrease of the Elliott- Yafet Ti 
in silicon structures with reduced spatial inversion sym- 
metry (gated, for example), since in such cases the inter- 
ference could be gradually removed. A recent experiment 
indeed finds a dramatic reduction of the spin relaxation 
time in Si/SiC>2 interfaces compared to the bulk |26| . 

Finally, we explore the anisotropy of Ti with respect 
to the valley orientation. Figure Bid) gives the ratio 
Tix/Tiz resolved for the intra and intervalley scatter- 
ing. Assuming that the anisotropy of Ti comes mainly 



from the spin mixing, that is, T Y ^ cx \J (b 2 )T,i(b 2 )T.j , 
we get Tif™/Tif™ « 2, as well as T^ ter /T/| ter « 
4— 2\/2 for intervalley /-processes (non-opposite valleys), 
and T*x ez /Tiz ter ~ 2 for ^-processes (opposite valleys). 
These estimates agree with our numerical results. The 
predicted anisotropy could be tested in strained silicon 
with lifted valley degeneracy. 

In summary, we explained the measured spin relax- 
ation in silicon performing realistic pseudopotential cal- 
culations of the Elliott- Yafet mechanism. We found the 
dominant scattering processes and their temperature de- 
pendence, predicted the valley anisotropy as well as Ti 
for hot electrons. 
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